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Abstract 

For the 2-brane Randall-Sundrum model, we calculate the bulk geometry for strong gravity, 
in the low matter density regime, for slowly varying matter sources. This is relevant for 
astrophysical or cosmological applications. The warped compactification means the radion 
can not be written as a homogeneous mode in the orbifold coordinate, and we introduce 
it by extending the coordinate patch approach of the linear theory to the non-linear case. 
The negative tension brane is taken to be in vacuum. For conformally invariant matter on 
the positive tension brane, we solve the bulk geometry as a derivative expansion, formally 
summing the 'Kaluza-Klein' contributions to all orders. For general matter we compute the 
Einstein equations to leading order, finding a scalar-tensor theory with oc ^/(l — \E r ), 

and geometrically interpret the radion. We comment that this radion scalar may become 
large in the context of strong gravity with low density matter. Equations of state allowing 
(p — 3P) to be negative, can exhibit behavior where the matter decreases the distance between 
the 2 branes, which we illustrate numerically for static star solutions using an incompressible 
fluid. For increasing stellar density, the branes become close before the upper mass limit, 
but after violation of the dominant energy condition. This raises the interesting question of 
whether astrophysically reasonable matter, and initial data, could cause branes to collide at 
low energy, such as in dynamical collapse. 
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1 Introduction 

Much progress has been made in understanding the long range gravitational response of branes 
at orbifold fixed planes, to localized matter. The Randall-Sundrum models, with one or two 
branes [1-3] simply have gravity and a cosmological constant in the bulk, making the linear problem 
very tractable. Linear calculations [4-7] show that in the one brane case, long range 4-dimensional 
gravity is recovered on the brane. Little is known of the general non-linear behavior [8,9]. There 
is no mass gap, and thus the non-linear problem is essentially a 5-dimensional one, even for long 
wavelength sources. Second order perturbation theory [10,11], and fully non-linear studies [12,13] 
are again consistent with recovering usual 4-dimcnsional strong gravity. In particular [12] shows 
that non-perturbative phenomena, such as the upper mass limit for static stars, extend smoothly 
from large to small objects, whose characteristic sizes are taken relative to the AdS length. 

For the two brane case, linear theory [4] found that the effective gravity is Brans-Dicke, for 
an observer on the positive tension brane, and a vacuum negative tension brane. For reasonable 
brane separations phenomenologically acceptable gravity is recovered. For observers on the neg- 
ative tension brane it was found the response was incompatible with observation for any brane 
separation, although stabilizing the distance between the branes does allow one to recover stan- 
dard 4-dimensional gravity [14-20]. In this paper we consider an observer on the positive tension 
brane, the negative tension brane to be in vacuum, and the orbifold radius to be unstabilized. 
This allows us to study the dynamics of the radion in strong gravity. Of course, our methods may 
be extended to include the stabilized case too. Although moduli have previously been assumed 
to be fixed at late times, a crucial feature of the recent cyclic Ekpyrotic scenario [21,22], is that 
the radius of the orbifold is not stabilized, the cyclic nature of this model being such that the 
separation always remains finite. 

The two brane strong gravity case appears more tractable than for one brane, as the linear the- 
ory allows a Kaluza-Klein style reduction of the propagator [4]. In Kaluza-Klcin compactifications, 
only the homogeneous zero modes are excited on long wavelengths, and the matter is thought to 
comprise fields over the whole internal space. However, for matter to be supported on the orbifold 
branes there must be modes excited which are not homogeneous in the extra dimension. Thus 
one cannot simply write down a non-linear ansatz for the metric, as is familiar from Kaluza-Klein 
compactifications. 

In previous work on the Horava-Witten compactification [23] , the low energy, long range effec- 
tive theory was constructed using consistent orbifold reductions [24-28]. The methods developed 
treated the homogeneous component of the metric as a background, and inhomogeneous pertur- 
bations about it as the contribution of the massive 'Kaluza-Klein' modes. The metric was then 
solved to leading order in a derivative expansion. Strong gravity was not discussed in these works 
except when considering inflation [28]. In order to use these constructions, all the zero modes of 
the orbifold must be homogeneous. One important result of the warped compactification, is that 
whilst the graviton modes are homogeneous, the radion zero mode is not [29]. This appears to 
be a general feature of warped models [30-32] , and one cannot directly use these Horava-Witten 
methods. 

The aim of this paper is to apply the orbifold reduction to the warped Randall-Sundrum model, 
and then examine the non-linear behavior of the radion. We begin in section 2 by illustrating the 
reduction proposed in [24-28] with a field theory example and discuss its application to strong 
gravity. In order to apply the method we must find a way to include the inhomogeneous radion 
zero mode. Instead of explicitly using an ansatz including the radion, we include it by 'deflecting' 
the brane relative to a non-linear extension of the 'Randall-Sundrum gauge' [4,5]. In sections 4 
and 5 we solve the bulk metric for conformally invariant, low energy, long wavelength matter on 
the positive tension brane. For conformally invariant matter, this coordinate system is Gaussian 
normal to the brane. The bulk geometry is constructed using a derivative expansion, which is 
formally summed to all orders. 

In section 6, we extend the linear analysis to the non-linear, for general low density, long 
wavelength matter, allowing the brane to become deflected relative to this coordinate system. This 
scalar deflection becomes the radion, and is treated non-linearly. This is crucial, as this field takes 
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values roughly of order the Newtonian potential, for static systems. Therefore, the case of strong 
gravity is exactly where the radion must be treated non-linearly. We calculate a local action for the 
effective induced Einstein equations, to leading order in the derivative expansion. The resulting 
theory is a scalar-tensor model, which reduces to Brans-Dicke in the linear approximation. This 
is the first non-linear, covariant derivation of the effective action. Considering the zero modes [33] 
does not tell one the conformal metric that matter couples to. Whilst cosmological solutions such 
as [15, 17, 34, 35] may derive radion dynamics, this only applies for a homogeneous radion, and 
thus is not a covariant derivation. Indeed the examples [15, 17] illustrate this point well. The 
form of the metric chosen, whilst coincidentally reproducing the correct effective action [33] , does 
not solve the correct linearized equations [29], and thus only applies to the case of cosmological 
symmetry. 

We discuss in detail the validity of the approximation for considering low density, long wave- 
length strong brane gravity. Terms that are quadratic in the energy density are neglected in 
the orbifold reduction method, consistent with the low energy approximation. The characteristic 
length scale of the matter must be large compared to the compactification scale for the derivative 
expansion to be valid. For strong gravity, the approximation then holds provided the 4-dimcnsional 
induced brane curvature invariants are small, compared to the compactification scales. Thus this 
method will not allow a global solution to a black hole geometry, or any other spacetimes with 
curvature singularities. However, it will apply to all other cases of strong gravity, such as static 
relativistic stars, an example being neutron stars, or dynamical non-linear systems, such as binary 
neutron star systems, or collapse of matter up to the point where curvatures become singular. 

Having pointed out that the radion may take large values for strong gravity configurations, we 
illustrate this using the example of relativistic stars. Previous work [36,37] has considered stars on 
branes using the projection formalism of [38], allowing the quadratic stress tensor corrections to be 
calculated by making ansatzes on the bulk geometry. However, these corrections are assumed to be 
negligible in our low energy density assumption, and it is the bulk geometric corrections which are 
relevant, and cannot be calculated in the projection approach. In [12] the full bulk solutions were 
numerically constructed for the one brane Randall-Sundrum case. Using the methods here, we are 
able to analytically construct the geometry for large stars in the 2-brane case. In section 9, using 
the leading order action, wc numerically consider the static relativistic star, for incompressible 
fluid matter. The linear theory shows that for perturbative stars with positive density the branes 
arc deflected apart. However, it indicates that if p — 3P becomes negative the opposite could occur. 
Now understanding the radion and bulk geometry non-linearly, we consider this, finding that as 
non-linear effects become important the branes do indeed become closer. Furthermore, the branes 
appear to meet before the upper mass limit is reached. This will occur for any brane separation. 
However, for phenomenologically acceptable separations [4], we find it does not occur before the 
dominant energy condition is violated. Neutron stars are believed to have polytropic equations 
of state that do not support negative p — 3P. It then remains a very interesting, and tractable 

1 + 1 problem to understand whether dynamical systems may cause branes to collide, for realistic 
matter and initial data. We then have the possibility that physical matter at low energies and 
curvatures, compared to compactification scales, might cause brane collisions requiring a Planck 
energy physics description. 

2 Orbifold Reduction: A Field Theory Example 

We illustrate the method of consistent orbifold reduction using a simple scalar field theory. The 
method was originally used in the context of Horava-Witten reductions [24-28] . We take the scalar 
equation, 

□ d * + (^*) 2 + aft + (d y vf = o (i) 

in a Minkowski (d + l)-dimensional bulk with signature —,+,+,... and coordinates y, with 
p. = 0, 1, ... (d — 1). We consider a finite range for y, choosing units such that < y < 1. The 
operator Dd is the Laplacian on rf-dimcnsional Minkowski space formed from x^. The analogy of 
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a Kaluza-Klein compactifications is to use periodic boundary conditions in y, identifying the field 
at y = and y = 1. For orbifold brane compactifications the boundary conditions to consider are 
Neumann, with the field gradient d y ^ = p(x) at y = and being zero gradient at y = 1. This 
will later correspond to the case of matter on the positive tension brane, and a vacuum negative 
tension brane. The zero modes of the system, with p(x) = are simply, 

tf(jc,y) = tf(jc), where, D d tf + (d^) 2 = (2) 

Having no y gradient this induces no p(x) on the boundary, and in brane language, is analogous 
to the zero mode ansatz discussed in [29]. Of course this solution also solves the Kaluza-Klein 
periodic boundary conditions, as these are identical to the Neumann ones if p(x) = 0. Then the 
full non-linear equations are solved in terms of a lower dimensional system. However, for non- 
trivial p(x) the solution cannot be independent of y, and the question is whether one can reduce 
the problem to a c?-dimcnsional one. 

If | p \« 1 we can linearize giving d d ^ + d 2 ^ = 0, whose solutions can obviously be found 
exactly. For slowly varying matter, we may use a derivative expansion in O d as, 

v) = + (4 + y \ y2 ) + w 

where wc implicitly assume that the c?-dimensional boundary conditions, relevant for the particular 
problem, are taken into account when evaluating the inverse Laplacian. It is these boundary 
conditions that specify the zero mode component in the solution 'P. However we now see that for 
large enough sources, of characteristic spacetime scale L and density p, the leading term will be 
large for pL 2 <~ 1. Then the non-linear terms in the original equation cannot be neglected. For 
low density matter, | p \« 1, we require L » 1 for such strong gravitational effects, and then 
we expect higher terms in the derivative expansion to become smaller. Note that if L < 1, the 
problem can be tackled simply using the linear theory, as | "J" |<< 1. 

For sources with L ~ 1 or smaller, nothing can be done in the non-linear regime. The problem 
is essentially a (d + l)-dimensional one. However, for L » 1, in the above derivative expansion, 
we see the leading term has no y dependence. One expects the remaining terms to be small, 
and thus we see might hope the non-linearity is not truly (d + l)-dimensional, but rather simply 
rf-dimensional. We now explicitly see how to realize this. 

2.1 Orbifold Reduction 

We now illustrate the orbifold reduction technique employed in [27, 28] , the issue being whether 
one can do better than perturbation theory for slowly varying sources. To proceed we observe that 
the leading term in the expansion (3), whilst being large, is independent of y. Therefore instead 
of linearizing the equation about ^f(x, y) — as above, we try separating a homogeneous piece of 
the solution from the inhomogeneous part, as, 

y(x,y) = H(x) + (f>{x,y) (4) 

where we choose that <fi(x, 0) = to define the splitting. The aim is to absorb the leading term 
of (3) into H(x). Then <f>(x,y) will consist of the remaining terms in the expansion (3), but these 
are all small. Note that H(x) is a zero mode solution only when p(x) — 0. The equation becomes, 

a d H(x) + (d^Hf + n d( j> + 2(d^H)(d^) + {d^f + d 2 y <p + (d v <pf = o (5) 

and assuming | <f) \« 1, can be linearized to eliminate quadratic terms in (f>, 

U d H(x) + (d^Hf + U d( j> + 2(d„H) + 8 2 <f> = (6) 

Removing such quadratic terms is the crucial step which will reduce the problem to a d-dimcnsional 
one. 
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2.2 'Strong Gravity' 

In order to characterize the magnitude of the terms we define, 



for any integer p > 0. We characterize, 



1 

a 



p(x) |~ L 2 p = $ 



(7) 



(8) 



which for large L allows $ ~ 0(1) even for small p. We now assume that we can perform the split 
so that, 



H(x) 



^- + h t + h 2 D d + ... 
Lid 



p(x) ~ 0(*) = 0(L 2 p) 



<f>{x, y) = {h(y) + f 2 (y)a d + . . . ] P {x) ~ O(p) (9) 
where hi are constants and fi(y) are only functions of y. Then we assess the terms above as, 



^(L 2 p) = 0(p) (j;(L 2 P )) 2 =0(<P P ) (i(L 2 p))(i(p)) = 0(p 2 ) O(p) O(^p) 

and now we see that in fact the cross term, (d fl H)(d^(j)) is of order 0(p 2 ), and can again be 
neglected. We are left with the equation, 



u d H + (d^Hf + d 2 y <p = -n dq 



(10) 



where terms on the left hand side are ~ O(p), 0(p<&) and on the right are order ~ 0(p/L 2 ). We 
also see that H and <j> have decoupled. The boundary condition for the problem is that, <f) = at 
y = and d y 4> = at y = 1. We may solve the system by taking some H(x) and then solving the 
linear problem, 



U d <t>(x, y) + d 2 0(x, y) = -D d [H{x)\ 



(11) 



where D d [H{x)\ = O d H(x) + (d^H) 2 and provides a homogeneous source term for <f>, of order 
O(p), even when H ~ 1. This equation can be solved exactly using a Greens function, and then 
derivative expanded as L >> 1 as, 



y-- 2 y 



£> d [ff(a;)] 



(12) 



Note that whilst we can solve this equation exactly, it is only useful if L » 1, as otherwise terms 
in the expansion will not decrease in magnitude. However for L < 1, $ << 1 in any case, and 
linear theory can be used. 

Now we have solved the (d+l)-dimensional problem as &(x,y) — H(x) + (j>{x,y), with the 
source p(x) = dy^ \ y =o, 



p(x) 



i + in d + o(n 2 ) 



D d [H{x) 



(13) 



and the non- linear (d + l)-dimensional problem is reduced to a d-dimensional one, which was the 
purpose of the exercise. Note also that we have not inverted Laplacians, and thus all issues of 



5 



Strong Brane Gravity and the Radion at Low Energies 



Toby Wiseman 



d-dimensional boundary conditions are implicit in the reduced equation. Specifying these, one 
may then invert the problem, 



D d [H{x)} = 



1 



n d + o{u\) 



p(x) 



= p(x) + KK corrections (14) 
where the zero modes are solutions for p{x) = 0. 

2.3 Regimes of Interest 

There are then two regimes of interest to us. Firstly the regime where $ ~ 1. In this case the sub- 
leading term in the derivative expansion Odp{x) ~ 0(p/L 2 ) is of order 0(p 2 ). Thus in this strong 
gravity regime, only the leading term in this derivative expansion is of relevance as we have not 
calculated the other 0(p 2 ) corrections due to linearizing in <j>. However the homogeneous split has 
allowed us to solve the non-linear problem up to 0(p 2 ) corrections. This is the case analogous to 
strong gravity on the brane for large objects, where effective corrections to the Einstein equations 
will be of order 0(l 2 K d p 2 ), with I the AdS length, n 2 the d-dimensional gravitational constant, 
and p the characteristic matter energy density. These are extremely small, and not of relevance 
in the strong field regime. 

The second case is where $ << 1, the weak gravity analogy. In this case we must contrast non- 
linear terms of order 0($>p) with terms in the derivative expansion, of order 0(p/ L 2p ). The former 
is a purely d-dimensional correction, whereas the latter depends on the size of the object compared 
to the fundamental length scale of the compactification. Note, in this example we chose units so 
that this fundamental scale was one. In the brane case, we are concerned with the ratio L/l, with I 
the AdS length. For large L, the non-linear terms in H{x), <~ 0($p), will be large compared to the 
derivative corrections. Simple linear theory would eliminate these non-linear terms in H(x), and 
therefore this method allows usual d-dimensional non-linearity to be automatically included. Of 
course the linear Greens function solutions tell one about all sizes of object and therefore contain 
vastly more information than the above method can yield. However, in cases where one is simply 
using the linear theory to derive the bulk metric in the long wavelength regime, it is far more 
powerful to use the techniques here. 

Note that for large enough p, one expects the derivative expansion term 0(p/L 2p ) to be of 
order 0{p 2 ), and then further terms no longer give meaningful corrections without calculating 
non-derivative 0{p 2 ) corrections too. The number of terms in the derivative expansion that are 
relevant before 0(p 2 ) is reached depends on the exact values of p, L. Certainly the first sub-leading 
term, p = 1, is always important if $ << 1, having magnitude 0(p 2 /$). 

3 Non-Linear Metric Decomposition 

In following sections we use the orbifold reduction to solve the non-linear field equations for an 
orbifold with matter on the positive tension brane, and a vacuum negative tension brane. The 
solution will relate the (d+l)-dimensional geometry to the usual d-dimensional non-linear Einstein 
equations. We take the (d + l)-dimensional bulk metric, 

I 2 

ds 2 = -^(g^(x, z)dx^dx v + dz 2 ) (15) 

with Greek indices taking values over the brane spacetime d-dimensions. If g^ v {x) is the Minkowski 
metric then this is simply AdS in Poincare coordinates, supported by a bulk cosmological constant 
K 2 d+1 K = —d{d— 1) / (2Z 2 ) , solving the bulk Einstein equations, 

G AB = (16) 
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where A, B are (d + l)-dimensional spacetime indices and g AB is the metric as in (15). In 
vacuum, positive and negative tension orbifold branes can be supported at z — z\ , z 2 respectively, 
with z\ < Z2- Observers are taken to reside on the positive tension Z 2 orbifold brane in addition 
to localized stress energy. We leave the negative tension brane in vacuum simply for convenience. 
However, all the methods outlined in this paper can be applied relaxing this condition. Note 
that the linear theory [4] already shows that an observer on the negative tension brane sees a 
phenomcnologically unacceptable gravity theory unless the orbifold is stabilized [18]. The orbifold 
Z 2 planes and localized matter I T I1V are treated in the thin wall approximation, the planes having 
tensions n 2 d+1 <y = ±2(d — I) /I. 

Following the methods of [24-28] , illustrated previously for the scalar field example in section 
2, we decompose g^ v into a homogeneous and inhomogeneous piece with respect to the coordinate 
z, as 

g^{x, z) = °g^ v (x) + h^ u (x, z) (17) 

There is obviously freedom in the above decomposition. To uniquely define it, we require that 
g^vix, z\) — °g^ v (x) for some constant z\, which will be the position of the positive tension brane 
if the matter has vanishing stress energy trace. Then the intrinsic metric on this surface Z\ will 
just be °g llv {x)l 2 1 z\. If the metric °g liv is Ricci flat, the (d + l)-dimensional Einstein equations 
are solved for vanishing [39]. These geometric deformations are then the gravitational zero 
modes of the orbifold. 

As discussed, we are not able to include the radion zero mode in the background, and have it 
remain homogeneous, due to the warped geometry [29]. Thus we only include the homogeneous 
gravitons, °g^ v {x). Instead we will include the radion by extending the brane deflection ideas 
of the linear theory. Thus we do not perturb the zz metric component as we will include any 
perturbative back-reaction of the radion in the non-linear deflection of the brane. 

The procedure we outline below is to find a suitable °g^ v {x) such that h^ v {x, z) remains small 
for a low density long wavelength matter perturbations on the orbifold brane, even when the 
intrinsic geometry is non-linear. Thus we aim to absorb all the non-linearity into °g^(x), which 
effectively shapes the induced geometry on the brane. One can think of h^ v as the contribution 
of the massive KK modes, which provide only small corrections to the induced geometry on the 
brane, but play the essential role of supporting the localized matter. 

As in the scalar field example in section 2, we define two dimensionless quantities. The first, p, 
characterizes the matter density or curvature scale compared to the AdS energy density. The sec- 
ond, L, compares the length scale associated with the matter, to the AdS length. For convenience 
we consider the parameter e = l/L 2 which is small in the relevant large object limit. Formally we 
take, 

P = l 2 \\'R^ a p\\ 



pe = l 4 \\'n>R^ a0 \\ (18) 

and then expect, 

peV~O{\\{l 2 >Uyi 2 >R» v a0 \\) (19) 

where p > and || || indicates the maximum absolute value of the tensor over all space, with I, the 
AdS length, introduced to make p, e dimensionless. The curvature tensors are formed from the 
induced brane metric, as is the Laplacian r O, and have indices arranged as they would appear in 
curvature invariants. In order to linearize the Einstein equations in the Kaluza-Klein contribution 
h^v, we will require that the matter be low density compared to AdS scales. In addition, in order 
to use a derivative expansion, we will require the object to have a large size or dynamical time 
compared to the AdS length. Thus we make the requirements that, 

p, e « 1 (20) 
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As mentioned previously, for p << 1 and e > 1, the problem can be solved simply using standard 
linear theory. An important point is that this restricts attention to non-singular geometries, as 
the curvatures must remain bounded, and our approximation will work only when they are well 
below the AdS curvature scale l/l 2 . For consistency we take the induced metric, t^„, to be of 
order, 

\\W 2p g^\\ ~0(^e*>) for p>l (21) 

Of course, when considering ||<7 M „||, the Ricci fiat zero modes must be taken into account. Thus 
we characterize the quantity, 

* = ||A&„||~0(f) (22) 

where we understand A<? M „ to be the difference of the metric from some Ricci flat zero mode 
background due to the presence of matter. We note that $ will be of order the Newtonian 
potential for static field configurations. 

As a technical note, we will also consider the magnitudes in equation (19) to hold, with the 
tensor indices arranged differently. This implicitly assumes that the metric is non-degenerate and 
non-singular. 

4 The Bulk Metric 

In this section, we solve the bulk geometry using a coordinate system which extends the Randall- 
Sundrum gauge to the non-linear case [4,5]. As in the linear theory, this coordinate system is 
such that constant z hyper-surfaces have scalar induced extrinsic curvature equal to d/l. We take 
a positive tension brane with localized matter at z = z\, and a vacuum negative tension brane 
at z = Z2, subject to the restrictions (20), and the further condition that zi,Z2 ~ 0(1). The 
geometry will be consistent only for positive tension brane matter with a vanishing stress energy 
trace. In the subsequent sections we remove this restriction, bending the brane relative to the 
surface z = z\ 1 again in analogy with the linear theory, although now these 'radion' deflections 
may be large, but slowly varying. 

Using the metric (15) and zero mode decomposition (17), the Einstein equations can be lin- 
earized in the Kaluza-Klein perturbation h^. As indicated in the scalar field example of section 
2, we will find that ~ ||^zfyui/|| ~ O(p) and therefore in linearizing we neglect terms of 

order 0(h 2 /I 2 ) ~ 0(p 2 /l 2 ) 1 consistent with the low curvature condition. Away from the branes, 
the linearized Einstein equations, with critical bulk cosmological constant, are, 

g zz n 2 d+1 Agii +1) = -r- ^- d * h + °( h2 n 2 ) 

G Zfl = \ (°V Q cU% - O V„0 Z JO + 0(h 2 /l 2 ) 

G^ - K^AgV+V = G^u - \(d 2 z h^ - °g^d 2 h) + ^A(d z h^ - °g^d z h) + 0(h 2 /l 2 ) (23) 

where R = °g^ v R^ v , and R^ u , G^ v arc the Ricci and Einstein curvature of the metric g^ v (x, z) on 
constant z hyper-surfaces. Indices are raised and lowered with respect to °g^ v (x), h = h a a and 
°V is the covariant derivative of this homogeneous mode metric. Thus, we are using "g^v in the 
same way one uses a background solution in usual linear perturbation theory. However, here, a g iiv 
is not a solution to these equations in the presence of brane sources. 
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As Wh^W ~ O(p), we may decompose the rf-dimensional curvature terms. The Ricci curvature 
term in G zz becomes 



R = °g^R^ = °g aPo Rap - °Dh + °V Q0 V^ Q(3 - °R a p hJ + 0(h 2 /l 2 ) 
= °R- °Dh + "V a0 V p h afj + 0( P 2 /l 2 ) 



(24) 



Note that the term °R a ph a P is of order 0(p 2 /l 2 ), as both l 2a R a [i and h a p are of order O(p). We 
now see that the linearized constraint equations, G zz — K^ +1 Ag^ +1 ^ = 0(p 2 /l 2 ), G Zfl — 0(p 2 /l 2 ) 
are satisfied if, 

°R(x) = 



h a a (x,z) = o V a h a ^x,z) = 



(25) 



This can be recognized as the usual Randall-Sundrum gauge condition. We have simply imposed 
this about a homogeneous background °g liV , which satisfies the condition °G a a = 0. The remaining 
Einstein equation becomes, 



- \z^d z (^f ) = \U V + Q(p 2 /l 2 ) 



(26) 



We emphasize that a G^ v and are independent of z. Note that ||°G pl ,||, 1| ~ 0(p/l 2 ), 

and °Dh^ u ~ 0(pe/l 2 ). For slowly varying matter, e << 1, and we construct the solution as a 
derivative expansion in "□. For our brane configuration, with matter only on the positive tension 
brane, we solve this differential equation using the ansatz, 



h^(x,z) =z*f(z, n) G flv (x) 



(27) 



which indeed satisfies the constraints (25) if °R = 0. The transverse condition is simply the 
contracted Bianchi identity, and the traceless condition holds as °R is zero. Then (26) reduces to 
the formal 'operator' equation, 



dll + -d z f 

z 



d 2 \ 



which can be solved exactly to give, 

2* (I)! 



(o D )f 3- ' * '21(2-1)! 
for arbitrary functions A(°D), B(°D), which can be expanded as, 

A(°n) = A a + A 1 a + ^-(°n) 2 + ... 



(°n)3iVd(\/ o nz) 



(28) 



(29) 



(30) 



and similarly for B, ensuring that the whole expression only contains positive powers of °D. The 
particular integral is taken as, 



/pi = < 



2 



2 



Z 2 



Z 2 



22(|-l)! 

-^--^n)i (n^^Uz) - iio g («D)j | (v^O. 



d = odd 



d = even 



(31) 



which again can be expressed as a derivative expansion in of Taylor series form. Then for long 
wavelength matter, defined through the condition (20), we expect leading terms to dominate the 
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series. By expanding the Bessel functions for v"Oz >> 1 in the particular integral solution (31), 
we see that the d-dimensional zero mode no longer well approximates the geometry, exactly as 
one predicts for small objects which should behave in a manifestly (d + l)-dimensional manner. 
If \AEl >> 1, ic. for small objects, e > 1 and $ = p/e « 1 so linear theory can be used. 
However, the case where it is z that is large, means one cannot study large non-linear objects in 
the one brane Randall- Sundrum case, as one must remove the second brane to distances greater 
than z > l/\AO. This essentially shows why the one brane case is really a (d + l)-dimensional 
problem. 

We must now fix the two functions of integration, A(°D), B(°E\), using the boundary condi- 
tions of the problem. At z = z 2 , the position of the negative tension vacuum brane, we require 
dzh^ix, z 2 ) — 0. The second boundary condition is simply the requirement that h fll/ (x,zi) = 
at z = Z\ by definition of the metric splitting (17). Direct evaluation of the expressions (29) and 
(31) for these conditions then determine A, B. 

Using this formal solution, represented by the operator C{z, °D) below, we may expand h^ v in 
powers of °D acting on "G^, the term (l 2o D) p °G M „ being of order 0(pe p /l 2 ), so that, 



V = C(z, °D) G^ = [h (z)l + h^zfU + o(°a 2 )] 



'G 



fJ,L> 



with, ho(z) — Bo — 



d-2 



A z a 



hi(z) 



,B z 2 



+ 2iim A ^ d+2 + Bi- 1 L - 2 z 2 + A l z d d = 3,d>5 



2(<2-2)- LJ U* 2(d-2)(d-4) 2(d+2)^ 

k \B a z 2 + ±A z 6 - i +lo f ~ 7 z 4 + |(logz)z 4 + B l - \ z 2 + A^ d = 4 



and, Aq 



1 



d(d - 2) z . 



d-2 ■ 



B« 



d-2 



Z2 



d-2 



(32) 



where Ai,Bi are trivially calculated from the boundary conditions but are not shown explicitly 
here for clarity. Whilst we have expanded to sub-leading order, O(l 4 ( D) G M i / ), as discussed in the 
scalar field example, in section 2, it is only the leading term that is relevant for strongly non-linear 
configurations with <I> <~ 0(1). Then the sub- leading correction to is of order 0(p 2 ), and one 
must also calculate the corrections from the linearization of the Einstein equations in d z h which 
enter at the same order. The sub- leading corrections are important in the case when $ < 1, and 
the field is weak, but one wishes to take e?-dimcnsional non-linear terms into account, such as when 
calculating post-Newtonian corrections. 



5 Conformally Invariant Matter 

For a positive tension orbifold brane at constant z = Z\, the induced metric is, 

p 



(33) 



as h^u = at z — z\ by construction. For an orbifold brane with tension cr, the Israel matching 
conditions [40] yield the localized stress energy 'T^ to be, 

~ K d+i <J 1 + K d+i Ir ^^ = ^ [K^v — 1 9tivK] z=Zi 



(34) 
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O(p) and the trace, h, vanishes. The terms proportional to 'g^ v cancel for the critical 
value of tension, a, leaving the localized matter, 

4+i T ^ = -- (d z C(z,°D)) \z= Xl °G^(x) (35) 

where we have substituted our solution for from equation (29). We see the trace of the stress 
energy tensor vanishes, as our construction has °R = so no trace can be supported on a brane at 
constant z. A vacuum orbifold brane with negative tension —a, is placed at z = Z2, where d z h^ v 
is zero. As the induced metric on the positive tension brane is simply given by equation (33), the 
induced Einstein tensor on the brane I G iiV = °G fll ,(x), and the Laplacian of the induced metric is 
J D = "Dz 2 /l 2 . Thus the effective Einstein equations for conformally invariant matter become, 

'G < x) = _^±il II 'T {x) 



—^ K 2 d) 'T^(x)+0(ep/l 2 ) 



with, «2 __-k2 and; I T a a = (36) 



2 1 

\d) - — y~ y^+i)' 

where we have formally summed all 'Kaluza-Klcin' contributions in the first line, and given the 
leading order approximation in the second. Indices are now raised and lowered with respect to 
the induced metric on the brane. To leading order in the derivative expansion this is simply d- 
dimensional Einstein gravity for conformal matter. This agrees with linear analysis, but of course 
is not restricted to linear matter configurations. 

In reference [41], the sub-leading correction has been written in the form of higher derivative 
terms in the action for the case of 2 branes stabilized by a bulk scalar. This was achieved by 
comparing such terms with the linear theory propagator. Here we have summed all terms in the 
derivative expansion, and explicitly used the non-linear theory. It would therefore be interesting 
to understand whether these may all be written in such a higher derivative form, as conjectured 
in this reference. 

The constant Bo in (32) was derived from the requirement that h^v vanish at z = z\. How- 
ever, we observes that in general it simply gives a small redefinition of the homogeneous metric 
component. Provided Bo is of order O(l), we might allow to be non-zero on the brane itself. 
However, the fact that 'G^ = a G^ u (x) exactly is a direct consequence of h^ u vanishing on the 
brane. If was not zero there, we would have to include this in calculating I G ilv . We now see 
the utility of this choice in simplifying the derivation of the Einstein equations. When considering 
non-conformal matter, this is the chief source of complication when discussing sub-leading terms 
in the derivative expansion, as in Appendix B. 

6 General Matter and the Leading Order Effective Action 

By identifying an extension of the Randall-Sundrum gauge, we have solved non-linear low energy 
gravity, providing 'T a a vanishes. The usual <i-dimensional response to conformal matter is re- 
covered on a brane at z = z\, with small higher derivative corrections, which we have formally 
calculated exactly. However we wish to allow the localized matter to have non-vanishing 'T a a . In 
the linearized theory one sees that placing the brane at non-constant z generates exactly such a 
trace. We therefore repeat this procedure non-linearly. We will find that the scalar field describing 
the deflection of the brane relative to the z coordinates, provides exactly the correct radion degree 
of freedom to 'replace' the trace of the matter stress energy, as happens in the linear theory. Note 
that the vacuum negative tension brane will remain at z — z^. Of course, using the same methods, 
one could extend the analysis to consider matter on the negative tension brane too. 
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For general matter we will only consider the leading term in the derivative expansion. Whilst 
we formally solve this to all orders in the previous section for conformally invariant matter, the 
non-conformal case is considerably more complicated. For clarity we only calculate using the 
leading order term, so the bulk metric is, 

I 2 

ds 2 = — [(°g^(x) + h^(x, z)) dx^dx" + dz 2 ] 



^g^(x) + (b - j^z 2 + A z d ^j °G^(x) + 0(°D°G)) dx»dx v + dz' 



(37) 



with Aq,Bo as in (32), although we note that this calculation could be extended to include the 
sub-leading terms in D°G, if required, as discussed in Appendix B. 

In order to calculate the matching conditions, we coordinate transform the metric to a new 
Gaussian normal system, z, adapted to the brane, which is located at z = z~\. This is in analogy 
with the linear theory, although we now consider large, but slowly varying brane deflections. 
Starting from the metric (37), we perform the coordinate transformation, 

z = ze^ s) , a;" = x» + (x, z) (38) 

where 4> is only a function of and plays the role of the d-dimensional scalar field parameterizing 
the deflection. We may choose the functions £ M to vanish on the brane for convenience. The 
technical derivation of the matching conditions is given in the Appendix A. Here, we simply 
outline the main results. The field <j> is treated non-linearly, although again its derivatives are 
taken to be small for large, low density objects. To leading order in the bulk derivative expansion, 
we find the induced metric on the brane is simply a conformal transformation of a g^ v , 

'9^{x) = ^e- 2 ^°g^(x) + 0{p) (39) 

where we note that x^ = x^ on the brane. The extrinsic curvature of the leading order metric 
(37) after coordinate transformation is then calculated at z = z\, and following on from that, the 
brane matching conditions. We eliminate °G MI , by tracing these, as in the linear theory, yielding 
the equation for the deflection, 

< nm + l_^v<^)) 2 = -2&=r$^r 1Ta <*w + ° {p2/l2) (40) 

which is, in fact, correct to all orders in the derivative expansion, ~ 0(e p p/l 2 ) for p > 0. Thus the 
radion receives no 'Kaluza-Klein' corrections as in the linear case [29]. By comparing the induced 
stress energy with the Einstein curvature of the above induced metric, eliminating °G^ V we may 
then derive the effective Einstein equations. Finally we are able to construct a local, covariant 
effective action for the leading order Einstein equations, in terms of a (i-dimensional radion scalar 
field, and the intrinsic metric on the brane. Then defining the following, as in Appendix A, 

m*)] = [- ir - ) , * d = - 2 — - 

A,M.x)] - CV/V, - 'g^D) <t>{x) - ('V^X'V^)) - ^-('W^x)) 2 ' 9fw (41) 

allows us to write the effective Einstein and radion equations (64) and (40), to order 0(max(^ 2/ D / T P i / , p 2 
as, 

(1 - n^(x)})'G^(x) = k 2 'T^ - (d- 2)n[cf>(x)]A fl „[<t>(x)] 

'4 = -A a a [<j>(x)} (42) 
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which result from the variation of the local action, 

S= 2^j ddX ~ (JV ^ )2 ) + / ^matterl'^] (43) 

where, 

m = 1 - fif>] = 1 - (le^)" 2 , and, = {^j^ (44) 

and —7} I T l _ lv 5 I g fl1 ' = <5£mattcr[ r S/ji/], with the matter Lagrangian having no dependence on the 
scalar (j>. The approximation of terms is fully discussed in the Appendix A. In addition, we discuss 
the calculation of sub-leading terms in the derivative expansion in Appendix B. 

We see that this is a general scalar-tensor theory for non-linear perturbations of <j>. For linear 
fluctuations we find, 

«M-£Kc!>"-i)+ <« 



= | [e 2y l l - l) , for rf = 4 (45) 

where y is the proper separation between the unperturbed branes. Thus for linear perturbations 
we find Brans-Dickc, with coupling as in [4]. Non-linearly the action is not Brans-Dicke. The 
same phenomenological gravitational constraint applies to the brane separation as in the linear 
theory [4], namely that y > 41. 

The form agrees with the non-linear zero mode action of [33], which was derived excluding 
matter and only including the orbifold zero modes. Strictly speaking, the orbifold zero modes 
do not show which conformal metric the matter couples to. Of course our treatment includes 
the graviton, matter and radion together, non-linearly. In [15, 17] the effective action was again 
derived through cosmological considerations, and again agrees, although, as discussed in [29,33], 
the radion ansatz used is incorrect, although does give the correct effective action for cosmological 
symmetry. 

There has been much work considering non-linear phenomena in scalar-tensor theory, with a 
view to observable tests [42]. Here we have the explicit form of the non- linear effective theory, 
and it will be interesting future work to explore its non-perturbative implications. One important 
point remains, namely the geometric interpretation of the scalar as the radion. If strong gravity 
phenomena, either static or dynamical, lead to a large, but order 0(1) <f> being generated, the 
branes could be forced together. This is discussed further in section 9. 



7 Geometric Interpretation of the Radion Scalar Field 

Having a consistent, non-linear realization of the low energy radion, allows us to formulate 
its 'meaning' geometrically, without reference to a particular background or coordinate system. 
Firstly, the positive tension brane is deflected to a position z = Z\e^ x ^ relative to the metric 
(15) of section 4, for general matter, the negative tension brane remaining at z = z^. The proper 
distance along the z axis between the two branes is, 

d(x) = l— = Zlog(4) - l<f>{x) (46) 

Jz ie 4> Z Z 1 

and thus the difference in proper distance due to matter is simply Ad(x) — —l(f>. In the linear 
theory one can define this radion as determining the proper distance between the branes, relative 
to the unperturbed background. Non-linearly, the coordinate system of section 4 is Gaussian 
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normal to the negative tension brane, and thus the z axis itself is a geodesic normal to this brane. 
Therefore, we find a geometric interpretation of the radion on the positive tension brane, at a point 
which relates its value to the proper length of the normal geodesic to the vacuum negative 
tension brane, that intersects the positive brane at x^. 

For non-conformal matter on the negative tension brane, so that this too must be deflected 
relative to the z coordinates, we again expect to find the radion to be related to the proper distance 
between the two branes along the z axis. More generally a line along the z axis can be phrased 
geometrically as a geodesic, normal to surfaces with constant induced scalar extrinsic curvature, 
of value d/l, related to the brane tension. 

8 Cosmology 

We now make contact with the previous result that some cosmologies can simply be constructed 
by considering moving branes in AdS [34,43-48]. This arises in our construction, through "g^ 
being either Minkowski or Milne spacetime, for such cosmological evolutions. This implies that 
"Gfu, vanishes, and with it, all the non-local correction terms. Thus, as one expects, cosmology is 
simply a result of the radion, and not the graviton zero modes or 'Kaluza-Klein' corrections. 

In order to see this, we write the induced FRW metric in conformal time, as a conformal 
transformation of Minkowski spacetime, 



However, we now observe from (39) that the metric °g^ u is similarly related to 1 g^ by a conformal 
transformation, 



to leading order in the derivative expansion. Therefore °g^ v must be a time dependent conformal 
transformation of Minkowski space. The remaining constraint is that °G flv is traceless, and there- 
fore the Ricci scalar curvature of °g llv must vanish. This restricts °g^ to simply be Minkowski 
space, or a Milne universe, and in both cases a G ilv vanishes, and the bulk is simply AdS. We have 
suppressed the fact that the metric above, (48), is only true to leading order in the derivative ex- 
pansion. However all the corrections to it are functions of °G M „, and so the leading order solution 
that °g llv is Minkowski spacetime, is true to all orders. 

Thus we have re-derived that fact that a subset of cosmologies considered in [34, 43-47] can 
simply be thought of as motion of a brane in pure AdS, which can be seen more generally for 
the AdS-Schwarzschild case, via a coordinate transformation [48,49]. Cosmologies based on AdS- 
Schwarzschild, when the black hole mass is not small, cannot be seen in this framework as they are 
not small deformations of the static Randall-Sundrum geometry. However, the machinery above 
could certainly be applied about an AdS-Schwarzschild vacuum solution with suitably modified 
brane tensions, and would remain largely unchanged, allowing all the generalized cosmologies to 
be included. 

We now briefly consider the validity of our approximations at late times. Then the brane scale 
factor is large, and therefore the positive tension brane is very near to the boundary of AdS, ie. 
z ~ 0, and thus <f> — ► — oo. So <j> ~ 0(1) is certainly not true. However, V</> = O(Lpl) ~ H ~ y/pl, 
with H the Hubble constant, which is still compatible with L 2 p ~ O(l). Thus, whilst <j> goes 
to negative infinity, the quantity z — ze^ simply moves more and more slowly to zero and the 
approximation still holds for late time cosmology. In fact we see that at very late times, when 
z — > 0, the effective scalar-tensor parameter becomes very large, indicating usual Einstein gravity 
is recovered. 




(47) 




(48) 
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9 Large Brane Deflection: Stars, Dynamics and Brane Col- 
lision 

In this section we consider whether branes can be forced together by low energy matter. In 
the previous sections, we have outlined the non-linear, geometric significance of the radion field. 
In addition, we have shown it to take values of order ~ O(l) for strong gravity configurations. 
Since phenomenological constraints have shown the branes may be minimally separated such that 
Zilz\ ~ e 4 , it may be possible that for very relativistic stars, or in dynamical collapse, the branes 
may actually touch, even though all curvatures remain low compared to the AdS length. It is 
instructive to consider the 4-dimensional linear theory for a static geometry, and perfect fluid 
matter, so that, 

'VV(x) = -f'T=|(p-3P) (49) 

Positive values of <j> will deflect the branes together, so we expect such behavior for 'T > 0. 
However, for usual matter 'T < 0, indicating that the branes in fact deflect apart for perturbative 
sources. We might expect this for non-linear sources too. Of course this intuition is based on the 
linear theory, and we must solve the full non- linear equations to confirm it. 

We numerically integrate the leading order equations from the action (43), for static spherical 
symmetry. Such configurations are relevant for relativistic stars, such as neutron stars. We do 
not present an exhaustive numerical study, but merely some illustrative examples, the solutions 
being intended to highlight the geometric nature of the scalar radion. We use a 4-dimensional 
spherically symmetric, static ansatz, 

ds 2 = -e 2T ^dt 2 +e 2R ^dr 2 +r 2 dn 2 2 



4> = 4>(r) (50) 

with dil 2 the area element on an S 2 . There have been several numerical studies of stars in Brans- 
Dickc and 'quadratic' models [50-54]. These studies all use polytropic equations of state, thought 
to be relevant for neutron stars. The effect we wish to illustrate is the dependence of the sign of 
the deflection on the sign of 'T = —p + 3P. Thus we choose the less realistic incompressible fluid 
matter, where 'T changes sign as non-linear effects become important. In standard GR, the core 
pressure diverges at the upper mass limit of the star, and thus *T becomes arbitrarily negative 
there. 

The incompressible fluid is taken to have density p = p . In the leading order action, the AdS 
length /, only appears in the usual lower dimensional Planck constant k 2 . Thus the solution of 
the above variables for a particular radius star, may be scaled in the usual manner to a star of 
any radius. For convenience, we choose 4-dimensional Planck units so that r is one unit at the 
stellar edge. In these units, the density corresponding to the upper mass limit in standard GR is 
po = 8/3, and is not expected to be very different in this scalar-tensor gravity [55]. Whilst this is 
not Brans-Dicke, the methods employed to numerically solve the stellar geometry are essentially 
identical. Initial data at r = is provided as 0o = 0(0), Pq — P{0), and the solution is assumed to 
be regular. The equations are Taylor expanded about r = to derive conditions for the remaining 
variables. These are then integrated using a fourth order Runge-Kutta algorithm. For a given p , 
one has a two variable shooting problem in (f>o, Pq. One must tune Pq in order to have P (1) = 0, 
and thus vanishing fluid pressure outside the star. In addition, 4>o must be varied to obtain 
</>(oo) = 0. In order to estimate this value, we simply integrate the equations to r = 10. However, 
in a more detailed study one should solve the vacuum exterior equations as in [50-54] and match 
the fields at the stellar edge. Note that for the scalar-tensor theory at hand, the scalar must take 
non-trivial values when matter is present, and thus one does not have to consider the issue of 
spontaneous scalarization as in [52-54]. 

A further point is the issue raised in [37] , that the extrinsic curvature of the brane is discontin- 
uous for an incompressible fluid. Thus we may not expect (19) to be well defined for p > 0. Whilst 
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1 2 3 4 r 5 

Figure 1: Plots of 4> against r for 3 values of core density, in units where the stellar edge is at r — 1. For 
low core density, the deflection is negative, the branes moving apart. For high density, near to the upper 
mass limit, p m ax — 2.67, the deflection becomes entirely positive, and the brane separation is decreased 
from the unperturbed case. 

this is true, it does not effect our leading order solution. It will be evident in sub-leading terms in 
the derivative expansion, which require °\3°G fll , to be evaluated, and appear problematic if °G^ V 
is discontinuous as is the case here. The physical resolution is to assume that on small scales, the 
fluid equation of state is softened to remove the problem. This will have negligible effect on the 
leading order solution. The only effect will be in the sub-leading terms, the contribution being 
suppressed by powers of l/L so f t , where L so ft is the length scale of the softening. 

In figure 1 we present three cases of star, plotting <f> against the radial coordinate r. The brane 
separation is chosen to be 41 in proper distance, compatible with the linear theory phenomenolog- 
ical requirements. However, we find the same generic behavior for all separations. 

We see that the lowest density configuration, with p = 1.0, does indeed behave as one expects 
from the linear theory, the deflection being negative in 0, and consequently the branes being 
further apart at the stellar core, than in the vacuum case. However, as expected, when the density 
is increased, —p + 'iP becomes positive in the interior and the deflection starts to become positive, 
as for p = 2.2. For p = 2.5, even nearer to the upper mass limit (p max ~ 2.67 from standard 
GR), the value of the scalar field is entirely positive. Note that for the branes to touch, <f) — 4 is 
necessary, so all these configurations are still far from this condition. 

Figure 2 shows the dependence of the core value of the scalar <f), on the density of the star p - 
We see the linear negative deflections, and then the turn-around to positive deflection for dense 
stars. The core radion value, <j>o, apparently diverges in the positive sense as the upper mass limit 
is reached, although more detailed numerical analysis would be required to explore this region. 
We mark on this plot the point where p — 3P changes sign. In addition the line p — P = is shown, 
being the point where the dominant energy condition is violated. One sees a tiny positive core 
deflection just before this, but nothing significant. In the same figure, the right-hand plot shows 
the core pressure against the density po, with the same lines marked. We see to gain a deflection 
of 0o — 1, one requires extremely large core pressures P ~ 40, compared to the core density 
po = 2.64. The data indicates that arbitrarily large <f)o can be reached very close to the upper 
mass limit, although the form of static matter required to support such a geometry is extremely 
unrealistic. 

In summary, whilst the scalar radion may become of order ~ O(l) in strong gravity, we find 
that for static stars, only extreme forms of matter can support configurations where the branes 
would meet. It is then a very interesting question whether the branes may collide in a dynamical 
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Figure 2: Left-hand frame; Plot of core radion value against stellar density, for incompressible fluid stars 
with fixed angular radius. The radion appears to diverge in the positive sense close to the upper mass 
limit. A radion value of <j>o ~ 4 would signal the branes touch. Whilst this will appear to happen very 
close to the mass limit, we note that (f>o has only just become positive when the dominant energy condition 
is violated. Thus, it appears for realistic matter, branes will not touch in a static, regular geometry. It 
is then interesting whether realistic matter may cause localized brane collisions in dynamical situations. 
Right-hand frame; Plot of core pressure against stellar density. The dominant energy condition is plotted, 
together with po — 3Po, whose sign reflects the direction of the brane deflection. 

process involving realistic matter. Such dynamics have been studied for Brans-Dicke [56-61] and it 
would be interesting to repeat these studies with this scalar-tensor model and realistic matter, to 
investigate whether there are initial configurations that will collapse to give positive <p large enough 
to cause a brane collision, in a region which is outside the formation of an apparent horizon. If such 
events could occur, it would indicate that the high energy string theory description of colliding 
branes, might be required to understand low energy, low curvature dynamics in astrophysical 
processes. 

10 Conclusion 

We have extended the linearized analysis of the Randall-Sundrum 2 brane compactification to 
allow the bulk, and induced geometry, to be calculated in the case of low energy strong gravity, 
with an unperturbed, vacuum negative tension brane. The induced geometry must have low 
curvature compared to the compactification curvature scale, and the matter energy density must 
be slowly varying on scales much larger than the compactification length. This is the situation 
relevant for non- linear gravitation in late time cosmology or astrophysics. 

The radion is inhomogeneous in the transverse coordinate for warped compactifications. It 
therefore can not be included as a homogeneous zero mode, as in usual orbifold reductions. We 
show how to extend the coordinate systems used in the linear analysis, to the non-linear, low 
energy case. Firstly we choose a coordinate system analogous to the linear 'Randall-Sundrum' 
gauge. This is Gaussian normal to branes with conformally invariant matter. We then solve 
the bulk geometry for conformally invariant matter on the positive tension brane, using a formal 
derivative expansion, which we sum to all orders. General matter is accommodated by introducing 
a non-linear 'deflection' of the positive tension brane with respect to this coordinate system. 
This deflection is governed by the radion, and is determined to all orders in the bulk derivative 
expansion. We are then able to understand the radion geometrically, as the proper distance 
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between the two orbifold branes, along a geodesic normal to the vacuum negative tension brane. 

Whilst we have applied the method to the 2-brane Randall-Sundrum geometry, the methods 
can be extended to any warped compactification where the radion is not a simple homogeneous 
mode. Note also, that we have considered a vacuum negative tension brane, allowing only matter 
on the positive brane. This is purely for convenience, and can be relaxed more generally using the 
same techniques. 

Having solved the bulk geometry, the effective action on the brane is calculated. To leading 
order, this yields a scalar-tensor theory with ui oc */(l — where \& is the scalar depending on 
the distance between the branes. We have noted that the radion scalar may take values of order 
~ O(l) in a strong gravity context, allowing the branes to touch or collide when gravitational 
fields become strong, and yet curvatures remain small compared to the compactification scale 
curvature. We have performed calculations for static spherical stars using incompressible fluid 
matter. For very low density stars, where the response remains linear, the branes are deflected 
apart. For densities where the gravitational response becomes strong, the brane separation may 
become reduced. We give evidence that the branes will touch before the stellar upper mass limit is 
reached, although the dominant energy condition is certainly violated. We find that, as indicated 
in the linear theory, p — 3P must become negative to get an appreciable effect. Whilst this does 
occur for incompressible fluid, it is thought not to for physically reasonable stellar material. Thus 
we find unrealistic matter is required to create touching branes in a static context. It raises the 
interesting question of whether realistic matter may allow branes to locally collide in a dynamic; 
context. If this were the case, high energy physics may be required to understand low energy 
astrophysical processes. 
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Appendix A: Non-linear Brane Matching Conditions 

We now give a detailed derivation of the effective Einstein equations for general matter, that was 
outlined in section 6. We start by considering the coordinate transformations on the bulk metric 
(37), 

z = ze^ s \ = x» +e(x,z) (51) 

Then preserving the Gaussian normal form of the metric requires that, 

d,e(x, z) = -ze+Wgr{x, z)B a e+& - gT{x, z)g ap {x, z)d a ?(x, z)d x ^(x, z) (52) 

where d a , d z are derivatives with respect to x^ and z. We see that there is the freedom to take £ M 
to vanish on the brane itself, as boundary data for this differential equation. The metric is then, 



I 2 

z l 



~ ms:) 9^(x, z)dx^dx v + N 2 dz 2 

g^(x, z) = g^(x, z) + g a ^(x, z)5^ a (x, z) + g af} (x, z)5^ a (x, z)d v £?{x, z) + z 2 e 24 ' ( ' x) d f _ t <j)(x)d u <j)(x) 

N 2 (x, z) = l + e- 2 ^g ap (x, z)(d z £ a (x, z))(d x ^(x, z)) (53) 

where we temporarily refrain from expanding the metric g^vix, z) in the coordinates x, z. 

From the linear theory [4, 5] we might expect (f> to be given by an expression similar to <fr ~ 
(l/U\)n 2 d T. Therefore <j> <~ <£> = O(pfe) and of course this cannot be assumed to be small. However, 
we do expect that gradients of <j) are small and use this to perform a gradient expansion of the 
transformed metric above, as was done for the bulk geometry in section 4. We formalize this with 
the assumption 

\\B a e*W\\ „ O(JL) , \\d a d^\\ ~ 0( P /l 2 ) (54) 

which we will see is consistent with the final result. Now considering the off diagonal equation 
(52), giving £ in terms of (f> 7 we see the quadratic term in £ is an order O(p) higher. We find to 
leading order, 

S z ^ = -ze+W°g l "*(x)d a e'K Si) [l + O(p)] (55) 

so that d z ^ <~ O(^j). We have decomposed g l _ lv (x,z) into a g^ u (x) and h l _ lv (x,z), and Taylor 
expanded in ^ = x^ — x^ . Repeating this for (53) we find, 

g^(x, z) = °g^(x) + [h^(x, ze+&>) + C*(x, z)B a °g^(x) + °g a(ll {x)d v) £, a {x, z) 

0(1) O(p) +z 2 (d^){d u e^)]+0{p 2 ) 

N 2 = l + e- 2 ^°g a/ s(x)d z ^(x, z)d x ^(x, z) + 0(p 2 ) 

0(1) 0(p) (56) 

Whilst we may Taylor expand h^ v {x + ^,ze^) in £, as £ is small ~ O(-^t), (f> is not small and 
may be of order one. However, we know the explicit z dependence of for each order in the 
derivative expansion and therefore we may simply evaluate /i M „ at z — ze^, rather than having to 
perform a Taylor expansion. Thus, 

V(*. ^ 0(2) ) = (b - ^"2^e 20(2) + A z d e d ^ °G^(x) + O(l 40 D a G^) (57) 
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As in the conformal matter case, /i M „ will not contribute to leading order in the induced metric, 
but will do in the localized stress energy. The induced metric on the brane at z = z~\ is simply, 

'g^{x) = ^e~ 2 ^g^(x, Z!) = ^e- 2 ^°g^(x) + 0{p) (58) 
z i z i 

where we use the fact that x a — x a on the brane. The O(p) terms, which involve h^ Ul are not 
relevant if one is only working to leading order in the derivative expansion, and their absence 
makes the following computations considerably simpler. 

We calculate the stress energy, K 2 d+1 'T^{x) = n^a'g^ + 2 [K^ v - K'g^]-^, as, 

K 2 d+1 'T^(x) = ~e-W*)(d xhlu/ faze+W)) \ s=Sl +21 A ^{x)] + 0( P 2 /l) 



A^ix)] = ('V/V„ - 'g^'D) <f>(x) - ('X7^(x))(<X7„4>(x)) t_l('V<t>(x)) 21 g^ v (59) 
where 'V is the covariant derivative of the induced brane metric 'g^(x) and, 

d z h^(x,ze^)) | f=fl = J^z ie 2 ^ + Aodzt'e^^ °G^(x) + O(/ 30 D°G) (60) 

As in the linear analysis, the trace of (59) determines <j>{x) as, 

V(*)(& K v {x, ze+&)) = ^e 2 ^0 z ( V(S) h^ix, ze+&)) + 0(p 2 ) = 0(p 2 ) (61) 

and hfu, is traceless with respect to "g^v This then determines the location of the brane in the z 
coordinate system by specifying that the radion <f> must obey, 

'U<j>{x) + d -^{'V4>{x)f = --J-^!$±±'r* a {x) + 0(p 2 /l 2 ) (62) 

This result is interesting as it is exact to all orders in the derivative expansion. From (58) the 
induced metric is simply a conformal transformation of the zero mode metric a 9p.v by e 2 ^ x \ to 
leading order, resulting in an induced Einstein tensor, 

'G^(x) = °G^(x) + (d- 2)A^ v [4>(x)} + 0(max(PD°G^,p 2 A 2 ) (63) 

where °\3°G flL , terms appear if one works to higher order in the derivative expansion including the 
terms O(p) in (58). All that remains is to solve (59) for °G^ U , and substitute into (63) to derive 
the induced Einstein equations on the brane, to leading order in the derivative expansion, 

'£„„(*) = ^ t - (d - ^"ftSL. -VM*)] + 0(max(i ! 'n«5'r„„, P 7i 2 )) 



0Wl)1 .(fi^)". 4 = ±^^ (64, 

which together with the equation (62) determining 4>{x), fully specifies the induced Einstein equa- 
tions, up to corrections O(l 2 °D G fJlU ) from sub-leading 'Kaluza-Klein' terms. Tracing these equa- 
tions, and using the equation for <f), one obtains, 

K 2 

'R = — *±±'T (65) 

which is true to all orders in the derivative expansion, the corrections being traceless. 

The sub-leading terms have a complicated form when written in the induced metric Laplacian, 
rather than "□, which is the Laplacian of °g ilv - However, of course these terms can be evaluated 
if required, and are briefly discussed in the following Appendix B. 
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Appendix B: Sub-leading Non-Local Terms 

We now consider the form that sub-leading non-local terms will take in the effective Einstein 
equations. Note that such corrections are only of interest when $ << 1, as otherwise they are of 
order 0(p 2 ), comparable with the other set of bulk corrections we have ignored. 

In the conformally invariant case we have formally summed the derivative expansion (32). This 
is considerably more complicated to do in the non-conformal case, the primary reason being that 
h^ v no longer vanishes on the brane, when it is not located at constant z as in the conformal case. 
Whilst the induced metric in the conformal case is simply 1 g^ v = (I 2 / ' z\)° g ^ , in the non-conformal 
case, 

= ^e- 2 ^ (<V(z) + [V0M1^ (X) ) + z 2 (d^)(d„e^)]) + 0( P 2 ) (66) 

which is presented above in equation (58) only to leading order as sub-leading terms were not 
considered. One can rewrite this as, 

Z ^e +2 ^'g^(x) = a g^{x) + [h^(x, f ie «*>) + ^("V^X^e^)] + 0(p 2 ) (67) 

and then calculate the Einstein curvature of both sides perturbatively using I g llv as the background 
metric on the left-hand side, and °g t _ iu on the right, giving, 

'G^(x) - (d - 2)A^[<P{x)\ = °G llv (x) - \°Uh^{x, ze+W) + 0( P 2 /l 2 ) (68) 

where the term z\ (°V p e^)(°V„e^) from the metric perturbation on the right hand side of (67) in 
fact contributes only to 0(p 2 /l 2 ). A subtlety in this calculation is that °V a h atJl (x, ze^) is not zero 
due to the implicit dependence. However, one finds that, 

°V a h aii (x, ze*) = (°V tt 0) zd- z h ali {x, ze*) ~ 0{(Lp)(^)) = 0(p 2 /l) (69) 

and so these terms can be ignored. Indeed, if $ ~ 1, then \\(j)\\ ~ 1 but the first sub-leading term is 
of order 0(p 2 ) and is not useful without calculating the other 0(p 2 ) corrections. If <f> << 1, then 
similarly ||0|| << 1. One finds that whilst it is important to consider the <fi dependence in the 
leading term of the derivative expansion of , which results in the recovery of the scalar-tensor 
behavior, for the sub-leading terms, <fr corrections are in fact negligible. Remembering we are 
considering the case ||0|| << 1, consider the two contributions, 

d>°G^ ~ (L 2 p)(p/l 2 ) 



</>°D°GV ~ (L 2 p)( — )(p/l 2 ) = 0( P 2 /l 2 ) (70) 

Whilst the first term, which could be found in the leading term of the expansion of h^ v is still 
relevant in the Einstein equations, the second term, found in the first sub-leading correction, is 
not as it is of order 0(p 2 /I 2 ). 

In order to write the effective Einstein equations as done in equation (64) , one must eliminate 
°G M1 / between (59) and (68). To leading order in the derivative expansion this is trivial as (68) is 
inverted to solve for °G^ V simply by rearrangement. However to sub- leading order this is a slightly 
more tricky procedure in general, as one must invert the derivative expansion. Of course this can 
be done if such terms are required, although we refrain from calculating them here. The final step 
is then to convert °G M „ to 'G^ and °D to 
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